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•^ I We examine the ability of LHC experiments to observe jets from squark pair production at 

•^ I the center-of-mass energy Ecm=10 TeV and 1 fb~^ of integrated luminosity. We point out the 

Q ' crucial influence of initial- and final-state radiation on the signal/background discriminating 

^— I . ability of different kinematic variables. The reliable measurements of missing transverse energy 

>• . and stransverse mass would play a key role in picking out the signal against the background. 
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1 Introduction 



Among many exciting challenges that the Large Hadron Collider (LHC) era 
brings most difficult is to build of the fundamental theory to describe the physics 
at arbitrary high energies. Supersymmetry (SUSY)|[T] is regarded as a widely 
favored candidate. 

SUSY with beautiful mathematical grace is good almost in all senses: the 
quadratically divergent contributions to the scalar sector are canceled, the light 
Higgs boson is predicted, the SUSY particle spectra contain the dark matter 
candidates. SUSY allows the unification of all known fundamental forces and, 
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presumable, allows explanation of baryon asymmetry in the universe, super- 
string theory naturally incorporates SUSY. 

To exclude some unwanted superpotential terms and provide the proton 
longevity the SUSY theories are usually supplied with discrete symmetry - 
R-parity conservation [1][J. As a consequence, the SUSY superpartners of the 
Standard Model (SM) particles are produced in pairs and decay to a SM particle 
and the lightest supersymmetric particle (LSP). The LSP is neutral, stable and 
weakly interacting and usually identified as the lightest neutralino or gravitino 
escaping the detector unseen. 

Unfortunately, the origin of SUSY breaking is not enough explored yet. The 
parametrization of soft SUSY breaking in supersymmetric Lagrangian density 
is now the basis of many models used in SUSY searches in collider experiments. 
Among them are the models with graviton, gauge and anomaly mediated soft 
SUSY breaking terms. Soft parameters can be determined from the LHC and 
ILC data if SUSY is discovered. 

The most popular SUSY breaking mechanism is realized in minimal super- 
gravity (mSUGRA)[5], where the SUSY breaking mediated from "hidden" to 
the "visible" sector by gravitational interactions. These models are strongly 
motivated, consistent with experimental data, and simple enough to explore in 
current collider experimentq3- In mSUGRA five parameters, four continuous 
and one discrete, 

are enough to calculate superpartner masses and mixings. Here mo is the com- 
mon mass for scalars, mi/2 is the common gaugino mass, ^40 is the common soft 
trilinear SUSY breaking parameter, sign[ii) is the sign of the Higgsino mass 
term, ton/? is the ratio of the vacuum expectation values of two Higgs doublets 
giving mass to the up and down type quarks. The parameters mo, mi/2, ^o 
are defined at the grand unification scale and tanP at the electroweak scale. 
Note that low-energy parametrization of the Minimal Supersymmetric Stan- 
dard Model contains up to 120 parameters arising mainly in the soft symmetry 
breaking terms. 

Intensive searches for manifestation of SUSY at the Tevatron were performed 
by CDF and DO collaborations (see e.g. Ref. [7] and refs therein). The LEP 



^The KK-parity conserving universal extra dimensional models and little Higgs[3] models with conserved 
T-parity as extensions of the Standard Model, predicts their own candidates for dark matter. 
^ About some problems in mSUGRA and its solutions see e.g. Ref. [6]. 



and HERA Collaborations also looked for SUSY in a variety of channels (see 
Refs. PB]). 



The ATLAS and CMS experiments have developed search strategies |IDl [H 
covering different SUSY breaking schemes and event topologies. 

At the hadron colliders the inclusive production of squarks and gluinos via 
strong interaction dominate over associated production of charginos and neu- 
tralinos and is one of the most promising discovery channels for SUSY. 

In this paper, we restrict ourselves to the study of the squark pair production 
in pp collisions at the center-of-mass energy Ecm = 10 TeV. This channel 
suffers from high backgrounds from QCD jet events. The requirement of high 
missing transverse energy values effectively rejects this background. We will also 
focus on the azimuthal angle between two hardest jets, the variables 777.^2 [12] 



and a [IB]. To eliminate multistep cascade decays of squarks we choose the 
parameter point of mSUGRA where squarks are lighter than the gluinos. In 
this case, the primary decay modes are q — > qx^. Thus, the squark pair will 
be presented in the interaction products as two quarks and two neutralinos. 
The final state exhibits missing energy and no mass peaks can be directly 
reconstructed. 

In addition to the QCD events, the main background to SUSY dijet signa- 
tures are electroweak processes such as VF + jets, Z + jets, tt + jets events. 
We find that the diboson production rate in the background is negligibly small 
for this study. 

The rest of this paper is organized as follows. The next section gives a brief 
description of signal and background events, the choice of the SUSY bench- 
mark point and some parameters of the generic LHC detector. Section 3 gives 
our strategy for the search for SUSY dijet events. We make use of different 
kinematic variables and examine signal and background events. In this Section 
we also make an estimate of the qq signal-to-background ratio. The impact 
of the calorimeter energy resolution, initial- and final-state radiation and jet 
reconstruction algorithm to the dijet variable distributions are studied. We end 
with the conclusions in Section 4. 



Signal and Standard Model background simulations 
for the generic LHC detector 



In order to study the ability of the LHC to observe the SUSY dijet events 
at EcM = 10 TeV we used the PYTHIA6.4||T1] event generator for signal and 
QCD background simulation. The ALPGEN2.12[T3] code was used for tt + nj, 



W + nj, Z + nj and diboson samples, where nj denotes n light jets. It was 
interfaced with PYTHIA for subsequent jet showering and hadronization. We 
include the following processes in our background analysis: 0-2 jets in tt, 0-3 
jets in Z and 0-4 jets in W events. Only in one diboson (ZZ, ZW, WW) event 
for 1 fb^^ of data set the sum of the transverse momenta of two hardest jets 
becomes greater than 500 GeV (see below). 

For this analysis, we used the mSUGRA scenario with the commonly used 
SPSla'[IB] benchmark point. The model parameters at this point are uiq = 
70 GeV, mi/2 = 250 GeV, Aq = -300 GeV, tan^ = 10, /i > 0. We used 
the SPhenofU] package for calculating the sparticle mass spectra, decay widths 
and branching ratios. The benchmark point used lead to a SUSY particle 
spectrum with m^i ~98 GeV, rriq ~550 GeV, rrig ~608 GeV. The total leading- 
order cross section of squark pair production at this point is ~4.65 pb if the 
transverse momentum of the outgoing partons in the hard scattering process is 
greater than 50 GeV. 

We use the central values of leading-order parton distribution function set 
from CTEQ6L1[IB1 and do not change the PYTHIA6.4 default choices for Q^ 



definition as well as factorization/renormalization scales. The initial- and final- 
state QCD and QED radiation (IFSR) and multiple interactions (MI) were 
enabled. 

The detector performance was simulated by using the publicly available PGS- 
4 [19] package written by J. Conway and modified by S. Mrenna for the generic 
LHC detector. The calorimeter granularity is set to (A0 x Ai]) = (0.10 x 0.10). 
Energy smearing in the hadronic calorimeter of the generic LHC detector is 
governed bjo 

^ = -^06 {Em GeV), 
where the stochastic term factor is a = 0.8 and the constant factor is 6 = 0.03. 



^We add the constant term to the PGS-4 simulation of energy smearing in the hadronic calorimeter. 



Jets were reconstructed down to \r]\ < 3 using the kr algorithm implemented 
in PGS-4. We chose D = 0.4 for the jet resolution parameter and required that 

1 2 

both leading jets carried a transverse momentum pj' > 50 GeV. 

In order to suppress backgrounds from the semileptonic SM and SUSY pro- 
cesses we select events without any isolated muon, electron, tau or photon with 
Pt > 20 GeV. We use the simphfied output from PGS-4, namely, a list of two 
most energetic jets. 

We simulated signal and background events at the rates corresponding to 
1 fb"^ of accumulated data. We note that the K-factor is ~1. 8-2.0 for QCD 
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Figure 1: Distributions of the E^^^^ (left) and if™'^^ (right) variables (see the text) for the 
signal and background events at Ecm = 10 TeV and integrated luminosity of 1 fb~^ and the 
SPSla' benchmark point. The sum of the transverse momenta of two hardest jets satisfies 
p^T > 500 GeV. 



Pt 



background , <1.2 for W+jets and Z+jets^, -1.1 for « + lj[ZI], and -0.89 
for ti + 2j|[22] events. The K-factor for the signal events at the benchmark point 
SPSla' calculated with using Prospino2.1[22] is about 1.5. In this analysis no 
K-factor was applied and we recognize that our results on signal significance 
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Table 1: The signal-to-background ratio for the qq production in pp interactions at Ecu 
10 TeV and 1 fb~^ of integrated luminosity estimated at different selection cuts for E^^'^^ . 



may be overestimated. 
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Figure 2: The normalized distribution of the difference of the azimutal angles for two hardest 
jets (p^ + Pt > 500 GeV) for signal events along with various Standard Model background 
sources at Eqm = 10 TeV and 1 fb"^ of integrated luminosity. 



3 Looking for squark pair production for one particular 
mSUGRA benchmark point 



Even in the early stages of its operation the LHC allows one to reach very 
large values of the jet transverse energy, the region which has never been studied 
before. 
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Figure 3: Expected distributions of the variables mT2 (left) and a (right) for the signal and 
background events at Ecm = 10 TeV and integrated luminosity of 1 fb~^ and SPSla' benchmark 
point. 



We start with comparing the missing transverse energy E™^^ of the whole 
event and H^'^^^ defined using only two leading jets in the event. In Fig. 1 we 
plot signal and background events binned in E^^^'"^ and H^''^^. The only selection 
criterion requires a hard cut of 500 GeV on the sum of the transverse momenta 
of two hardest jets. A comparison of the left and right panels of Fig. 1 shows 
that the effectiveness of background restrictions based on the Ej^^^'^ and iif^ *"'^* 
cuts would significantly differ towards high values of these variables. In Table 1 
we present the signal-to-background ratio, S/B, for different selection cuts for 
Ej^^^^. It would be possible to achieve the signal significancq^ S/y/S + B of 
about 10.5 with the accumulated data of 1 fb"^ using only the E^'^^^ cut. With 
only the cuts on H™^'^ it is impossible to get signal significance even close to 3. 

Besides being of physics origin, events with the missing transverse energy 
also have other sources. Mismeasurements of the energy of jets, incomplete 



'* We follow the tradition but can say that if no robust signature in the mass peak is observable, the use of 
this signal significance is not so compelling. 
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Figure 4: Expected i?™*** distributions for signal events at the SPSla' benchmark point (left) 
and QCD background events with and without initial- and final-state radiation, for different 
constant terms b in the hadron calorimeter resolution expression, and with the jet cone al- 
gorithm (right). LHC operation at Ecm = 10 TeV and integrated luminosity of 1 fb~^ is 
assumed. 



coverage of the calorimeters, large electronic noise, etc. would lead to a large 
^miss sig]^ai_ por this reason, in the early stage of the LHC running it may 
not be possible to use E^"^^^ as a signal/background discrimination variable. In 
what follows we do not implement any Ej^^^^ requirement. 

The dijet angular distribution is a useful measurement to probe the SM 
processes as well as new physics manifestation. In Fig. 2 we plot the normalized 
dijet angular distribution, 1/N dN/dA(f), where A0=0i — 02 for two leading jets. 
The signal event distribution is nearly fiat in contrast to the background event 
distributions. 

An important quantity which in principle allows observing squarks and deter- 
mining them masses and masses of their invisible decay products is the mrofl^j. 
Inspired by the transverse mass ttit in the W —> lu decay this variable is 
sometimes called the stransverse mass. For the decay of two massive identical 
invisible particles mT2 is defined as 

mT2{m)^ , min {max[m^(ri, ^^,^iv), m|(pj, ^|,^iv)]}. 
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Table 2: The signal-to-background ratio for the qq production in pp interactions at Ecu 
10 TeV and 1 fb~^ of integrated luminosity estimated at different selection cuts for mT2- 



Here p]p' are the sum of the momenta of the visible decay products of a parent 
particles, ^in is the trial mass parameter, namely, the mass of the LSP and 
mT{Pv,Pi, m,) = m^ + mj + 2{EyTE,T - PvtPzt), E,^v)t = \JpI{v)t + ^\v)- The 
minimization is taken over all possible missing energy jg^*** splittings. 

In Fig. 3 (left panel) we demonstrate the mT2 distribution for events with the 
sum of the transverse momenta of two hardest jets greater than 500 GeV. For 
this calculations we assume the massless LSP. The S/B for different mT2 cuts 
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Figure 5: Same as in Fig. 4 but for the difference of the azimuthal angles of two hardest jets. 

is presented in Table 2. The signal significance of 8.2 is achievable using the 
771^2 cut alone and ~10 with the combination of 771^2 and Ej^^^^. The tandem 
of H^'^^'^ and 171x2 does not lead to improvement of the signal significance. 
We find that 0.25(0.4)fb~^ of an integrated luminosity when the LHC runs at 
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Figure 6: Same as in Fig. 4 but for 772^2 distribution. 



^CM=10 TeV would provide observation of the qq signal with 5 standard devi- 
ations of statistical uncertainty using the E^^^^ {171x2) distribution. Further, we 
estimate that a factor of ~2.3 less integrated luminosity is required at 14 TeV, 
and a factor of about 2.5 more at 7 TeV, to achieve the same signal sensitivity 
as at 10 TeV. 

Recently, the new measurable a^^ has been proposed which can be used 
as a signal/background discriminating variable in squark pair production. This 
analysis does not explicitly require the E'^ *** cut. The variable a is defined as 
the ratio of the second hardest jet transverse momentum to the invariant mass 
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In Fig. 3 (right panel) we show the a distributions for the signal and back- 
ground events. As one can see, the clean SUSY signal at larger a values is not 
viable due to a bigger tail of the QCD background distribution. The a > 0.6 
signal is 45 times weaker than the SM background. 
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Figure 7: Same as in Fig. 4 but for a distribution. 



To get further insight, in Figs. 4-7 we show what happens with the Ej^^^'\ A0, 
171x2 and a distributions of the signal and QCD background events when IFSR 
are switched off in the simulated events. The impact of the constant term b in 
the calorimeter energy resolution formula and jet definition algorithm is also 
demonstrated in these figures. The cyan histogram in the figures corresponds 
to the cone jet definition algorithm with the cone radius R = 0.7. The signal 
event distributions are robust to all these influences. The presence of IFSR 
significantly increases the QCD background for all these distributions at large 
values of the E™'^'^, 111x2 and a and small values of AcJ). We can see that the 
missing transverse energy for the QCD background at large values of E^'^^'^ is 
less susceptible to the IFSR. The impact of the MI on the 111x2 distribution 
of the signal events is demonstrated in Fig. 8. One can see that the upper 
edge of these distributions indicates to the mass of the parent squark. The 
experimental resolution leads to the smearing of the correct mass point. 

Finally, we remark that the behavior of the a distribution of the QCD back- 
ground without IFSR in Fig. 7 looks surprisingly like the ones obtained in 
Refs. HSl 
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Figure 8: Expected 771^2 distributions for signal events with initial- and final-state radiation 
and multiple interactions and without them when the sum of the momenta of two hardest jets 
is greater than 500 GeV. The data correspond to 1 fb~^ of integrated luminosity at Ecm = 10 
TeV. 

4 Conclusions 



We have demonstrated that squark pair production in pp collisions at the 
center-of-mass energy of Ecm = 10 TeV is a promising channel for discovery of 
SUSY at the early stages of LHC running. The our attention was focused on 
the distinctive characteristics of background subtractions in qq production for 
the mSUGRA benchmark point SPSla' where the squarks are light. 

With 1 fb^^ of integrated luminosity, strong evidence for a qq can appear with 
hard cuts on the missing transverse energy and stransverse mass distributions. 
In the mSUGRA model with a benchmark point SPSla' it is possible to reach 
a signal significance 5a of statistical uncertainty with 0.4 fb^^ of accumulated 
data when LHC operates at Ecm = 10 TeV. 

We find that initial- and final- state radiation forms a crucial background for 
the a 113] distribution of the qq signal. It should be hard to pick out the signal 
against the background in SUSY dijet events using the measured a distribution 



even with a tandem of other measurables. The variable a not can provide an 
obvious advantage for discovery of qq as was stated previously. 

The careful validation of the whole richness of new physics at the LHC will 
be possible with improvement of our understanding of the known physics. 
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